Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Hallmarks of transcriptional intratumour heterogeneity across a thousand tumours.

Each tumour contains diverse cellular states that underlie intratumour heterogeneity (ITH), a central challenge of cancer therapeutics1. Dozens of recent studies have begun to describe ITH by single-cell RNA sequencing, but each study typically profiled only a small number of tumours and provided a narrow view of transcriptional ITH2. Here we curate, annotate and integrate the data from 77 different studies to reveal the patterns of transcriptional ITH across 1,163 tumour samples covering 24 tumour types. Among the malignant cells, we identify 41 consensus meta-programs, each consisting of dozens of genes that are coordinately upregulated in subpopulations of cells within many tumours. The meta-programs cover diverse cellular processes including both generic (for example, cell cycle and stress) and lineage-specific patterns that we map into 11 hallmarks of transcriptional ITH. Most meta-programs of carcinoma cells are similar to those identified in non-malignant epithelial cells, suggesting that a large fraction of malignant ITH programs are variable even before oncogenesis, reflecting the biology of their cell of origin. We further extended the meta-program analysis to six common non-malignant cell types and utilize these to map cell-cell interactions within the tumour microenvironment. In summary, we have assembled a comprehensive pan-cancer single-cell RNA-sequencing dataset, which is available through the Curated Cancer Cell Atlas website, and leveraged this dataset to carry out a systematic characterization of transcriptional ITH.

An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma.

Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Promiscuous antibodies characterised by their physico-chemical properties: From sequence to structure and back.

Human B cells produce antibodies, which bind to their cognate antigen based on distinct molecular properties of the antibody CDR loop. We have analysed a set of 10 antibodies showing a clear difference in their binding properties to a panel of antigens, resulting in two subsets of antibodies with a distinct binding phenotype. We call the observed binding multiplicity 'promiscuous' and selected physico-chemical CDRH3 characteristics and conformational preferences may characterise these promiscuous antibodies. To classify CDRH3 physico-chemical properties playing a role in their binding properties, we used statistical analyses of the sequences annotated by Kidera factors. To characterise structure-function requirements for antigen binding multiplicity we employed Molecular Modelling and Monte Carlo based coarse-grained simulations. The ability to predict the molecular causes of promiscuous, multi-binding behaviour would greatly improve the efficiency of the therapeutic antibody discovery process.

Significant Differences in Physicochemical Properties of Human Immunoglobulin Kappa and Lambda CDR3 Regions.

Antibody variable regions are composed of a heavy and a light chain, and in humans, there are two light chain isotypes: kappa and lambda. Despite their importance in receptor editing, the light chain is often overlooked in the antibody literature, with the focus being on the heavy chain complementarity-determining region (CDR)-H3 region. In this paper, we set out to investigate the physicochemical and structural differences between human kappa and lambda light chain CDR regions. We constructed a dataset containing over 29,000 light chain variable region sequences from IgM-transcribing, newly formed B cells isolated from human bone marrow and peripheral blood. We also used a published human naïve dataset to investigate the CDR-H3 properties of heavy chains paired with kappa and lambda light chains and probed the Protein Data Bank to investigate the structural differences between kappa and lambda antibody CDR regions. We found that kappa and lambda light chains have very different CDR physicochemical and structural properties, whereas the heavy chains with which they are paired do not differ significantly. We also observed that the mean CDR3 N nucleotide addition in the kappa, lambda, and heavy chain gene rearrangements are correlated within donors but can differ between donors. This indicates that terminal deoxynucleotidyl transferase may work with differing efficiencies between different people but the same efficiency in the different classes of immunoglobulin chain within one person. We have observed large differences in the physicochemical and structural properties of kappa and lambda light chain CDR regions. This may reflect different roles in the humoral immune response.